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Warp propagation in astrophysical discs 


Chris Nixon and Andrew King 


Abstract Astrophysical discs are often warped, that is, their orbital planes change 
with radius. This occurs whenever there is a non-axisymmetric force acting on the 
disc, for example the Lense-Thirring precession induced by a misaligned spinning 
black hole, or the gravitational pull of a misaligned companion. Such misalignments 
appear to be generic in astrophysics. The wide range of systems that can harbour 
warped discs - protostars. X-ray binaries, tidal disruption events, quasars and others 
- allows for a rich variety in the disc’s response. Here we review the basic physics 
of warped discs and its implications. 


Big whorls have little whorls, which feed on their velocity. And little whorls have 
lesser whorls, and so on to viscosity. 

- Adaptation of The Siphonaptera by Lewis Richardson 


Chris Nixon 

JILA, University of Colorado & NIST, Boulder CO 80309-0440, USA 
Einstein Fellow 

e-mail: chris . nixon@ jila . Colorado . edu 
web: http://jila.Colorado.edu/chrisnixon/ 

Andrew King 

Department of Physics and Astronomy, University of Leicester, Leicester LEI 7RH, UK 
e-mail: ark@astro . le . ac . uk 

web: http://www2.le.ac.uk/departments/physics/people/andrewking/ark 


1 



2 


Nixon & King 


1 Introduction 

In general gas possesses angular momentum and orbits any central massive object 
such as a protostar or a black hole rather than falling directly on to it. The test par¬ 
ticle orbits are not closed ellipses but form rosettes, as perturbations such as nearby 
stars or gas self-gravity mean that the potential is not perfectly Keplerian. These 
orbits intersect and so the gas shocks, causing dissipation of orbital energy through 
heating and radiation. However angular momentum must be conserved during this 
process. Therefore the gas settles into the orbit of lowest energy for the given angu¬ 
lar momentum - a circle. In most cases the gas does not all have the same specific 
angular momentum and forms a disc rather than a ring. More generally gas has angu¬ 
lar momentum with a spread of directions and settles into a warped disc, where the 
disc plane changes with radius. An initially planar disc may become warped through 
a variety of different effects. If the disc is misaligned with respect to a component 
of the potential these include Lense-Thirring precession from a spinning black hole 
(Bardeen & Petterson, 1975), gravitational torques from a companion (Papaloizou 
& Terquem, 1995) or torques from a misaligned magnetic field (Lai, 1999). An ini¬ 
tially aligned disc may also be unstable to warping through processes such as tides 
(Lubow, 1992) or radiation from the central object (Pringle, 1996). 

A planar (unwarped) disc (Pringle & Rees 1972; Pringle 1981) evolves under 
the action of a turbulent viscosity. For simplicity this is often parameterised in the 
Shakura-Sunyaev form where the viscosity is written as v = ac^H (Shakura & 
Sunyaev, 1973). The most likely process driving viscosity is the magneto-rotational 
instability (MRI) suggested by Balbus & Hawley (1991). This instability requires a 
weakly magnetised disc whose angular velocity decreases outwards, as is the case 
in Keplerian and near-Keplerian discs. The cartoon picture of the instability has 
two parcels of gas, at different radii, joined by a magnetic field loop. The inner 
parcel rotates faster, so the loop becomes stretched. The magnetic field acts to slow 
the inner parcel and speed up the outer parcel, thus transferring angular momentum 
from the inner parcel to the outer parcel. The net effect is that the inner parcel moves 
inwards through the disc and the outer parcel moves outwards. In reality this process 
drives turbulence in the disc, which then ultimately transports angular momentum. 

The Shakura-Sunyaev viscosity parameterisation is a simple but powerful and 
intuitively appealing way of characterising the angular momentum transport in a 
planar disc. As described in Pringle (1981) any turbulent eddies must have a length- 
scale smaller than the disc thickness H. Moreover any supersonic turbulent motions 
in the disc shock and rapidly dissipate, giving a maximum signal velocity of ap¬ 
proximately the sound speed Cj. So the maximal viscosity expected in a disc is 
Vmax = c^H, and we can write the viscosity with the Shakura-Sunyaev a parameter 
as V = acsH, but it is important to realise that a need not be constant in position or 
time. For fully ionised discs, observations suggest that a « 0.1 — 0.4 (King et al., 
2007). Numerical simulations of the MRI without net vertical field do not yet re¬ 
produce this result, instead finding a « a few x 10^^ (e.g. Simon et al., 2012). In 
other systems where the gas is neutral or only partially ionised, observations allow 
the MRI-driven a to be lower. For example, observations of protostellar discs find 
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a ~ 0.01 (e.g. Hartmann et al., 1998), but in dead zones this can be much lower 
(~ 10^“*, e.g. Fleming & Stone 2003; Simon et al. 2011; Gressel et al. 2012; Martin 
& Lubow 2014). 

In a planar disc, where the gas is essentially hydrostatic in the vertical direction, 
the evolution is usually described by a one-dimensional viscosity which commu¬ 
nicates angular momentum radially. However, in a warped disc the communication 
of angular momentum is three-dimensional. Early investigations (e.g. Bardeen & 
Petterson 1975; Petterson 1977; Petterson 1978; Hatchett et al. 1981) simply as¬ 
sumed that the viscosity could be described by the same v in all directions, and the 
evolution equations they derived did not conserve angular momentum. 

The first self-consistent investigation was provided by Papaloizou & Pringle 
(1983). They give two derivations of the equations of motion for a warped disc. 
The first (their Section 2) is close to the previous derivations in the literature, but 
gives the correct form of the internal torques required to conserve angular momen¬ 
tum. The second derivation (their Section 3) is the first to take the internal fluid 
dynamics into account. This derivation is the first to discuss the radial pressure gra¬ 
dients induced by the warp that oscillate at the orbital frequency, creating a reso¬ 
nance. Therefore the torque corresponding to the vertical shear (or more precisely 
the radial communication of the component of angular momentum parallel to the 
local orbital plane) depends inversely^ on the viscosity parameter, a, as the reso¬ 
nant forcing drives shearing motions which are damped by a (see Section 4.1 of 
Lodato & Pringle, 2007, for a discussion). An implication of Papaloizou & Pringle 
(1983) is that simple approaches to the dynamics of warped discs can miss subtle but 
important effects. But the complex nature of warped discs makes both simple and 
elaborate approaches necessary - for physical insight and to get the physics right. 
We shall see that the simple approach can capture most, but not all, of the relevant 
physics. 

The work of Papaloizou & Pringle (1983) formally requires H/R < a <C 1, and 
most restrictively, that the disc tilt (U^/RQ in their notation, where t/j is the per¬ 
turbed vertical velocity) he ^ H/R. This last restriction requires that the warp be 
so small as to be unobservable. So to make progress, Pringle (1992) derived an 
equation of motion using only conservation laws, and not the internal hydrodynam¬ 
ics, but retaining full generality in terms of the disc tilt. It is this approach that we 
will follow in deriving the evolution equations below (Section 2.1.1). These equa¬ 
tions are formally valid (in that they conserve angular momentum) for arbitrary disc 
warps and viscosities (Vi,V 2 ). However, this approach offers no insight into deter¬ 
mining the values of the viscosities or how they depend on any system parameters 
such as warp amplitude. In this sense its treatment of the internal fluid dynamics is 
“naive” (Papaloizou & Pringle, 1983). 

For this reason, Ogilvie (1999) started from fully three-dimensional fluid equa¬ 
tions (compressible, with a locally isotropic [‘Navier-Stokes’] viscosity), and de¬ 
rived a full evolution equation for a warped disc with these properties. This deriva¬ 
tion confirmed the equations derived by Pringle, with two important differences. 


* See the definition of the quantity A at the end of page 1189 of Papaloizou & Pringle (1983). 
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First, the Pringle equations did not include a torque that tends to cause rings to pre- 
cess when tilted with respect to their neighbours, and second, the torque coefficients 
are determined as a function of warp amplitude and a. The extra torque between 
rings is not required for angular momentum conservation, so it does not appear in 
the derivation in Pringle (1992) and the torque coefficients are determined by the 
local fluid dynamics between rings. Papaloizou & Pringle (1983) derived the torque 
coefficients for the azimuthal and vertical shear terms in the linear approximation, 
obtaining the well-known result that ai = a & ai = 1/(2a). Ogilvie (1999) ex¬ 
tends this to the fully nonlinear regime, determining the coefficients for arbitrary 
warp amplitudes. This work marked a great advance in the understanding of warped 
discs. Ogilvie (2000) extends these equations to include the effects of viscous dissi¬ 
pation and radiative transport. 

In cases where a < H/R the evolution of a warped disc is not controlled by 
viscosity. Pressure forces dominate and can propagate warping disturbances through 
distances > R as waves. This occurs because the a damping is too slow to damp 
the wave locally (on scales of order H) allowing it to propagate. A second, more 
subtle, condition must also be satisfied for wave propagation to be efficient - the 
disc must be close to Keplerian, i.e. |l — < H/R (see e.g. Section 5.1 of 

Wijers & Pringle, 1999). This is required to make the forcing resonant in a warp 
(Keplerian implies the orbital, vertical and epicyclic frequencies are the same, Q = 
Q^ = k). The first dedicated study of the behaviour of warped discs in this wave-like 
regime is Papaloizou & Lin (1995). This paper showed that these waves, driven by 
pressure gradients, propagate at approximately half the local sound speed. Further, 
the propagation behaviour becomes diffusive when a small viscosity a ^ H/R is 
included. Lubow & Pringle (1993) and Korycansky & Pringle (1995) investigated 
the propagation of a more general class of waves in discs. As discussed by Pringle 
(1999) it is possible to recover the warp wave propagation velocity = Cs/2 (for 
a = 0) from this approach. 

In this review we discuss the two types of warp propagation, through waves and 
diffusion. We discuss the evolution equations, their interpretation and derivation. We 
describe the viscosity, and in particular the relation between the small scale turbulent 
a viscosity and its role in shaping the effective viscosities which control the dynam¬ 
ics of warped discs. Finally we discuss some major results and some outstanding 
problems in understanding this complex and subtle accretion disc behaviour. 


2 Warp propagation 


In a planar accretion disc there is only azimuthal shear acting in the disc driving a 
turbulent viscosity which transports angular momentum. In a warped disc there is a 
second ‘vertical’ type of shear and so the rate of orbital shear is not simply parallel 
to the local disc normal /, but is given by (e.g. Ogilvie & Latter, 2013b) 


S = 


ds dQ dl 

R—=R l + RQ — 

dR dR dR 


( 1 ) 
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where s {R,t) = Q[R)l {R,t) is the local orbital angular velocity. We can think of 
the evolution of a warped disc as controlled by two torques, communicating angular 
momentum along and normal to the local orbital plane respectively. 

The usual planar disc viscosity is driven by azimuthal shear and is assumed to 
result in a turbulent a viscosity acting against the shear and transporting angular 
momentum radially. The second viscosity associated with the warp is induced by 
radial pressure gradients driven by the misalignment of neighbouring rings (see 
Fig. 1). In a near-Keplerian disc the forcing frequency (epicyclic) and the orbital 
frequency resonate to produce a significantly enhanced torque, whose response is 
controlled by a. This torque is driven by pressure and results in the launch of a 
pressure wave. When a > H/R this is damped locally and the evolution is diffusive, 
but when a < H/R the wave can propagate large distances. This leads to two types 
of warp propagation in discs; diffusive and wave-like. These are the subject of the 
next two sections. 

The convenient coordinates for describing a warped disc are a hybrid of cylindri¬ 
cal polars and Euler angles. A local annulus (of width ^ H) is described by cylin¬ 
drical polars but each ring of the disc is tilted in three dimensional space 

described by the (spherical) Euler angles j3 {R,t) & Y{R,t), which correspond to the 
local disc tilt and twist respectively. So we define the disc unit tilt vector as 

/ = (cos ysin p , sin ysin j3, cos j3). ( 2 ) 

The local angular momentum vector for the disc is L = ER^Ql, where E is the disc 
surface density. The (Cartesian) surface of the disc is 

X {R,(l)) = R (cos 0 sin 7 -b sin (j) cos /cos j3, sin 0 sin ycos j3 — cos (j) cos 7 , — sin 0 sin j3) , 

(3) 

where the azimuthal angle (j) is zero at the descending node and increases in the 
direction of the flow. Eurther quantities such as the elements of surface area then 
follow as described in Pringle (1996). Vector calculus in these warped disc coordi¬ 
nates is comprehensively summarised in Section 2 of Ogilvie (1999). We next derive 
the equations of motion for warped accretion discs. 


2.1 Diffusion 

Here we consider diffusive discs with ////? < a < 1. We derive the equations of mo¬ 
tion in a simple way, allowing significant insight into the problem. This derivation is 
naive to the internal fluid dynamics in a warped disc (Papaloizou & Pringle, 1983), 
and so we compare and contrast with the results of more complete derivations that 
are too involved for our purposes here (e.g. Ogilvie, 1999). 
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Fig. 1 This figure (cf Fig. 10 of Lodato & Pringle 2007) illustrates the radial pressure gradient 
indueed by a warp. The top and bottom panel show cross-sections of the same two neighbouring 
rings of gas, but at different azimuths (j). The shaded regions indicate the higher pressure around the 
local midplane, and the arrows show the resultant pressure gradient when the rings are misaligned. 
The azimuthal angle around a ring is measured in the direction of the flow from the descending 
node where (j) = 0. The tilted rings cross at the nodes, and so at these points are in aligned contact as 
usual. At all other azimuths the ring midplanes do not fully line up, causing a region of overpressure 
above or below the midplane. Each gas parcel feels an oscillating pressure gradient as it orbits in 
the warp. 


2.1.1 Evolution equations 

The derivation here is a three-dimensional version of the standard planar accretion 
disc equations (e.g. Pringle 1981; Frank et al. 2002), so we briefly recall the planar 
case first. 

We assume the disc is planar and that all quantities can be azimuthally averaged 
(requiring that the radial velocity is much smaller than the orbital velocity, Vs ^ V^) 
so all quantities are functions of radius R and time t only. The disc has surface 
density E {R,t), radial velocity Vr and angular velocity Q {R,t). We consider 
an annulus of gas between R and R + AR. The mass in this annulus is 2nRARE. 
Conservation of mass relates the rate of change of this mass to the net flow of mass 
into or out of the annulus 
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[InRARL) = {ItiRVrL)^ - {InRVRL)^^^;,. (4) 

Rearranging and taking the limit as AR —> 0 gives the continuity equation (mass 
conservation) for a disc 

Similarly we can derive the equation expressing angular momentum conservation. 
The angular momentum of an annulus is IkRARLR^Q. and the rate of change of 
angular momentum is the net flux of angular momentum plus the net torques; 

^ {inRARLR^Q) = {inRLR^ Q.Vr) - {inRZR^ ^IVr) (6) 
+G(R + AR)-G(R) . 


Again by rearranging and taking the limit AR —0 we have the angular momentum 
equation 


l(i:R^n) + l^(Ri:vRR^n) 


1 dG 
2%r'M' 


(7) 


where G{R,t) is the internal disc torque resulting from the disc viscosity. The vis¬ 
cous force is proportional to the rate of shearing 


F 


dQ 

= 27zRH^iR—, 


( 8 ) 


where Rdf2 /dR is the rate of shear and /r is the dynamic viscosity, related to the 
kinematic viscosity hy jj. = pv, where p is the density. Using p ^ E /// the force is 


F 


=InRvER- 


df 2 

T_ 

dR 


(9) 


Note that this force acts in the azimuthal direction, and thus when taking the cross 
product with the radial vector to generate the torque we see that the torque acts in 
the correct direction (i.e. on the z-component of angular momentum - we are only 
considering a planar disc so far). So we can write the internal viscous torque as 


G = InRvEREl'R 


( 10 ) 


and the equation expressing angular momentum conservation as 


dt 


{ER?Q.) 




(REVrR^Q) 


iA 

R'dR 


(yER^Q.') . 


( 11 ) 


Multiplying Equation 5 by R^El and then subtracting from Equation 11 we rear¬ 
range for Vr to get 


RE-^{R^n) ■ 


(12) 
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Substituting this into (5) gives 

dt RdR[ 


which for the Keplerian potential, Q = GM/R^, is 


dE 3 d 

~ rJr 


(13) 


(14) 


This is a diffusion equation, so if we know the viscosity v we know how a distri¬ 
bution E {R,t) will evolve. So far we have said nothing about the disc viscosity or 
what drives it - we will return to this in detail in Section 2.3. 

Assuming the disc is in vertical (z-direction) hydrostatic balance we derive the 
vertical structure by equating the relevant components of the pressure force and 
gravity 


I dP _ GM z 
p dz 


(15) 


where the last factor is the geometrical sine term. If the disc is thin we have |z| <C R, 
so 


1 dP _ GMz 
p dz R^ 


(16) 


Now if we assume an equation of state, e.g. isothermal P = c^p, then 


2 <3lnp GMz 


(17) 


and we can integrate to get 


where we have introduced the disc scale-height P[ = c^/Q.. From this we see that a 
disc in vertical hydrostatic equilibrium has H/R = Cs/V^, so for thin discs the flow 
is supersonic. 

We can now attempt a similar calculation for a warped disc, where the inter¬ 
nal torques act in three dimensions. This calculation draws heavily on Papaloizou 
& Pringle (1983) and Pringle (1992). The local angular momentum density is 
L{R,t) = ER^£2l{R,t), where l{R,t) is a unit vector in the direction of the ring 
angular momentum. The mass conservation equation is the same as for the flat disc 
(5). Angular momentum conservation is expressed as 
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^ (inRZR^QARl) = {ItzRZR^QVrI) (19) 

-{IkRZR^QVrI)^^^^ 

+GiR + AR)-G{R) , 


where G is the three-dimensional internal torque. Again rearranging and taking the 
limit AR 0 we have 

+ = ,20) 

The viscous torque has two obvious components. The (/?, 0) stress contributes a 
torque acting in the direction of / (cf 10 ) 

Gi=2nRviERn'Rl, (21) 


where Vi is the azimuthal shear viscosity. For two neighbouring rings with 1 and 
/ + Al, the {R,z) stress acts to communicate Al between the rings, so the torque acts 
in the direction dl/dR. The {R,z) torque is then 

G2=27zR^V2ER^Q^, (22) 


where V 2 is the vertical shear viscosity and the factor of a half comes from integrat¬ 
ing cos^ (j) (cf Fig. 1) across the ring (Papaloizou & Pringle, 1983). There is also a 
third component (Papaloizou & Pringle, 1983; Ogilvie, 1999) which we could write 
as 

G 3 = 27tRV2,ER^Ql X . (23) 

But this defines an effective viscosity V 3 - this is not an appealing notation as this 
term does not lead to diffusive behaviour. Instead this torque causes a ring to precess 
if it is inclined with respect to its neighbours. This produces dispersive wave-like 
propagation of the warp (Ogilvie, 1999). 

We combine (21), (22) & (23) to give the internal torque (cf. Equation 55 of 
Ogilvie & Latter, 2013b) 


G = 2nRER^Q 


/Q''' 

\ 1 dl 

dl 

i'"' 

)(+2V3^ + v,(x 



Now, putting this into (20) we have 


liER^ni) 


r'^ 


(EVrR^QI) 


LA. 

rJr 

LA. 

LA. 


(viER^Q'l) 

(v3Z«=£>;x|1), 


(24) 


(25) 
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which is identical to that found in Pringle (1992) except for the last term. 

Pringle (1992) showed that this equation can be combined with (5) to eliminate 
Vr. This is done by taking the dot product of / with (25) and then subtracting R^Q x 
Equation (5) and using I ■ dl/dR = 0 and d/dR{l- dl/dR) —0^1- dH/dR^ = 
— \dl/dR\^. This gives (Pringle, 1992, Equation 2.3) 


Vr 




\v2l.R^Q. 



RE^iR^Q) 


(26) 


We note that the extra torque derived by Ogilvie (1999) does not have a component 
in the direction of / (cf. Equation 123 of Ogilvie 1999). This shows its precessional 
nature: it causes neighbouring rings to change their planes, but does not drive any 
radial flux of mass or angular momentum. Einally, substituting Vr into (25), we get 
the evolution equation for the disc angular momentum vector L 


dt 


lA 

r'^ 


Ud/dR) ViER?(^-a' 


+ 


+ 


1 d ri 


E{d/dR){R^Q) 

dl 


^ -V2R\L\ ^ 

RdR 2 ' ' dR 


1 d \\\v2R^Vl\dl/dR\^ 
^^1 {d/dR) (R^n) 



\ d f T dl 
RdR\ ^ dR 


(27) 


We see that this equation evolves both the disc shape l{R,t) and surface density 
E (Rf) through the angular momentum vector L = ER^Ql. This equation is 
almost exactly that derived by Pringle (1992), with the addition of the last term de¬ 
rived by Ogilvie (1999). The last term can also be found in the linear hydrodynamic 
analysis of Papaloizou & Pringle (1983) where the warp diffusion coefficient is 
complex (see also Equation 2.1 of Kumar & Pringle 1985).^ This term leads to pre¬ 
cession in the presence of a warp (Ogilvie, 1999), but its coefficient, V 3 , is smaller 
than V 2 by a factor « a. Therefore in time-dependent problems this term is often 
neglected (e.g. Section 3.4 of Lodato & Pringle, 2007). However, in the case of an 
inviscid non-Keplerian disc, this is the only non-zero internal torque (Section 7.2 
of Ogilvie, 1999). 

If the rotation law and the torque coefficients (Vi,V 2 ,V 3 ) are specified in terms 
of the disc quantities, equation (27) can be solved numerically, conserving angular 


^ See the definition of the quantity A at the end of page 1189 of Papaloizou & Pringle (1983). In 
the linear equations it is often convenient to assume the unit tilt vector I = ( 4 , 4 , 4 ) has 4 ~ 1 and 
therefore adopt complex equations for the disc tilt where e.g. W = lx + ily In this case a diffusion 
coeffieient with non-zero real and imaginary parts has eomponents in the direction of both dl/dR 
and I X dl/dR. 
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momentum to machine precision by the method described in Pringle (1992). This 
has been used in many papers to explore the evolution of warped discs (e.g. Pringle 
1997; Wijers & Pringle 1999; Lodato & Pringle 2006; Nixon & King 2012). 

The derivation of this equation was somewhat simplistic, as it did not take the in¬ 
ternal fluid dynamics into account (Papaloizou & Pringle, 1983). So to check its va¬ 
lidity Ogilvie (1999) started from the full three-dimensional fluid-dynamical equa¬ 
tions and derived the equations of motion (Equations 121 & 122 in Ogilvie 1999). 
Ogilvie (1999) was able to confirm the equations derived by Pringle (1992), with 
two important refinements. First, the Pringle (1992) equations were missing the V 3 
term which makes neighbouring annuli precess if they are tilted with respect to each 
other (included above), and second, the torque coefficients are uniquely determined 
by the warp amplitude and a. 

For comparison we give the evolution equations derived by Ogilvie, written in 
our notation but retaining the Qi coefficient form. The mass conservation equation 
is 

dH \ d , 

^ + = o (28) 

and angular momentum conservation is 




-t- 


'r'm 


Qi-E^R^Ql 


c? 3 dl 
Q2T.—R^Q — 
^ Q dR 


QsE^R^Qlx-^ 
^ n dR 


(29) 


where we have used the relation ^ = Ec^/£2^ for the azimuthally averaged second 
vertical moment of the density.^ We should compare these equations with (5) & 
(25). The two sets of equations become identical given three relations between the 
torque coefficients (Vi,V 2 ,V 3 ) and the effective viscosity coefficients (2i, 02 , 23 ) 
(Fodato & Price, 2010; Nixon & King, 2012), i.e. 

Vi oc ai (a, Irl) = («, Irl) = -^2i («,|V^|), (30) 

V 2 °-a 2 (a,|r|) = 222(a,|r|) (31) 

and 

V3°=«3(«,|r|) = 23(a,lrl) , (32) 


^ Strictly speaking the rhs of this relation is missing a factor of order unity dependent on the warp 
amplitude | yr|, as the forcing in a warp restricts the hydrostatic balance assumed in its derivation. 
This affects the viscosity coefficients Qi,Q2,Ql significantly for large | \j /\, and so in this case we 
must use the form of Ogilvie (2000) taking this effect into account (Ogilvie, private communica¬ 
tion). 
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where |y/| = /? \dl/dR\ is the warp amplitude and the first relation is given in sim¬ 
plified form for the case of Keplerian rotation. 

The coefficients Qi, Q 2 & Qi depend on the warp amplitude |v^| = R\dl/dR\, 
a and the orbital shear (see Figures 3, 4 & 5 of Ogilvie 1999). In various limits 
their approximate values can be obtained analytically (e.g. Section 7.3 of Ogilvie 
1999) but in general need to be computed numerically. For a complete generalisa¬ 
tion of the Shakura-Sunyaev theory of planar discs to diffusive warped discs, see 
Ogilvie (2000) where the effects of viscous dissipation and radiation transport are 
also included in the viscosity coefficients. 

The equations derived in this section, and those by Ogilvie (1999, 2000), hold 
when a > H/R (damping the otherwise uncontrolled resonant response to the 
epicyclic forcing) and when H/R^ 1. However, the main uncertainty in modelling 
a ‘real’ astrophysical disc in this manner is the nature of the turbulence driving an¬ 
gular momentum transport. The approach adopted above is supported by strong evi¬ 
dence (analytical, numerical and observational) but it is certainly not beyond doubt. 
We comment on this in more detail in Section 2.3 below. 


2.2 Waves 

The evolution equations for a wave-like disc were first derived by Papaloizou & 
Lin (1995), Demianski & Ivanov (1997) and Lubow & Ogilvie (2000). All of these 
papers linearise the fluid equations, and there is very little work on the nonlinear 
behaviour of wave-like discs. Ogilvie (2006) extends the equations into the weakly 
nonlinear regime, but this has yet to be followed up in any detail. There is good 
reason to believe that the nonlinear dynamics may well display many interesting 
features (Gammie et al., 2000; Ogilvie & Latter, 2013a,b), but there remains no fully 
nonlinear theory for warp propagation in this case. However, the linear dynamics of 
wave-like discs is reasonably well understood and we discuss this here. 


2.2.1 Evolution equations 

The equations of Lubow & Ogilvie (2000) follow a similar notation to that used in 
Section 2.1.1 and so we shall use that here also. The linearised equations of motion 
for a wave-like disc are one for the conservation of angular momentum 


2 dl 1 dG 
ZR^Q— = -—+T 
dt R dR 


(33) 


and one for the evolution of the internal torque 


dt 


-h 




(.2 37 

IxG + aQG^ZR^Q^—. 

4 dR 


20 . 


(34) 
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The external torque T can arise from a variety of effects, but is zero if the disc 
is precisely Keplerian and there are no external effects. The second term on the 
Ihs of (34) is also zero when the disc is Keplerian {£2 = Q,= k). The third term 
on the Ihs describes damping of propagating waves, which occurs on a timescale 
^ 1/ {a£2) - we shall refer to this as a damping. If we compare (33) to (20) we can 
see that these equations make the assumption that E{R) is independent of time and 
the radial velocity is zero - i.e. there is no diffusion of mass or advection of angular 
momentum in these equations. 

In the simplest case of an inviscid {a = 0) Keplerian disc the equations become 


2 dl 1 dG 
dt R dR 

(35) 

— =ER^Q — —, 
dt 4 dR' 

(36) 


and 


which can be combined into a single equation by eliminating G 

1 d 


dH 

dt^ ~ dR 


4 dR 


(37) 


Equation 5 of Ogilvie (2006) shows that when EH is independent of R, this can 
be recast into a classical wave equation for the disc tilt. This shows that the wave 
speed is Vw = Cs/2 (Papaloizou & Lin, 1995). The wave speed can be evaluated in a 
similar way from equations A53-A56 of Lubow & Ogilvie (2000) to give the more 
general result that jAE. 

With the equations in this reduced form (35 & 36) it is easy to see where the terms 
come from. Equation (35) is a statement of angular momentum conservation in the 
case where E is independent of time and so = 0 (cf. Equation 20). Equation (36) 
expresses how the internal torque changes because of the pressure applied by a warp, 
which might be written as dGjdt = RH.pV^.R^Q.dl/dR, which with p = E/H 
gives (36). This implies that each ring of the disc responds as if hit by a pressure 
wave with velocity Vw. With such a simple approach to the internal torque we cannot 
recover the wave speed, but this allows physical insight into the disc response to 
a warp. By including a small a damping of the wave propagation, and any non- 
Keplerian terms due to the potential, we get the full evolution equations derived 
rigorously by Lubow & Ogilvie (2000). 

Taking the (artificial) limit that a is large in these equations (34) becomes 

T c? dl 

aQG = ER^Q — —. (38) 

4 dR 


Substituting this into (33) gives 


j dl \ d 
ER^Q— = - — 
dt RdR 


dl 


ER^Q 


4an dR 


(39) 
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Now we can see that in the linear regime with V 2 = cji jla this torque is exactly 
that derived to communicate the disc tilt in the diffusive case. This is not by chance; 
the physics of the two terms is precisely the same in this approximation. It is likely 
therefore that for specific (short) timescales one might use the linearised wave equa¬ 
tions with strong damping to model a linear diffusive disc (see e.g. Facchini et al., 
2013), but that this approach is ultimately missing much of the physics in warped 
discs (mass diffusion, advection of angular momentum and internal precession). 
Good agreement between the approaches may be found in cases where these effects 
are minimised by extra physics in the problem, for example where tidal torques 
inhibit mass flow. 

As a warp wave propagates in the disc it can lead to either a local or global 
bending of the disc. For a single wave the criterion for this is effectively governed 
by the wavelength (A) of the disturbance. For ^ A < the disc responds by 
warping, but when X > R the wave is unable to bend the disc which instead responds 
by tilting as a whole. This behaviour was observed in the simulations of Nixon & 
Pringle (2010) where a warp was induced at the outer edge of a disc by a passing 
perturber. This warp had a wavelength of approximately one third of the outer disc 
radius, leading to warped outer regions. As the warp wave propagates inwards the 
inner disc {R < A) tilts as a whole while remaining planar. 

In wave-dominated discs there is also the possibility that a disc might respond 
to a disturbance by globally precessing. In the paragraph above this clearly occurs 
when A > Rom- However, it can also occur when the wave communication across 
the entire disc is short compared to any precession time of the disc. In this case 
the disc is able to share the precession generated at each radius across the entire 
disc and respond as a cohesive whole. This has been observed in many simulations 
(e.g. Larwood & Papaloizou, 1997; Fragile et al., 2007; Fragner & Nelson, 2010), 
although there appears to be some ambiguity in the literature as to the exact criteria 
for the onset of this behaviour. It is likely that it requires waves to damp non-locally 
(effectively a <^H/R) and that the wave travel time across the disc (~ /?/cs) should 
be shorter than any precession induced in the disc. 

Nixon & King (2013) pointed out that when only a small piece of the disc is 
simulated these conditions can be artificially met, facilitating repeated global pre¬ 
cession. This occurs as the propagating wave reflects off the outer boundary, which 
in reality should be orders of magnitude further out. This allows unphysical global 
precession. Instead the wave should leave behind a steady disc shape (e.g. Lubow 
et al., 2002). In other cases, discs may be kept radially narrow by extra physics such 
as tidal truncation, allowing physical global precession (e.g. Lodato & Facchini, 
2013). Care must also be taken when simulating warped discs, as the propagation 
of waves depends sensitively on the ratio of a and H/R, these parameters must 
be accurately modelled for the system of interest. Often thick discs are employed 
in simulations to increase vertical resolution, without observational or theoretical 
motivation. Similarly the turbulent viscosity may not be modelled correctly if nu¬ 
merical affects (e.g. low resolution) lead to strong numerical dissipation, or if the 
effective a in simulations which resolve turbulence is smaller than that implied by 
observations (King et al., 2007). 
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The main untouched area for wave-like discs is a thorough understanding of their 
nonlinear evolution - Ogilvie (1999) discusses this for the diffusive case. In the 
wave-like case Ogilvie (2006) extends the equations of Lubow & Ogilvie (2000) 
to the weakly nonlinear regime, but the predictions made are yet to be tested by 
a hydrodynamics code. Also, large warps may be hydrodynamically unstable for 
small a. Gammie et al. (2000) predict that such warps are subject to a parametric 
instability which leads to dissipation and enhanced wave damping. More recently 
Ogilvie & Latter (2013a,b) have developed a local model for warped discs and have 
used this to explore the onset of turbulence in this hydrodynamic instability. Such 
instabilities have not appeared so far in global simulations, probably because of 
insufficient resolution and artificial dissipation effects. 


2.3 Viscosity 

We have seen in the previous sections that the usual approach to implementing a 
viscosity for warped discs is to adopt the a parameterisation, where the local stress 
is proportional to the pressure. This means that each component of the shear (hori¬ 
zontal and vertical) is damped by viscous dissipation at the same average rates. This 
is the origin of the term ‘isotropic viscosity’, but this does not imply that the effec¬ 
tive viscosities (torque coefficients Vi, V 2 & V 3 ) are isotropic. Instead, the internal 
structure of a warp means that these torque coefficients take very different values 
that depend on the disc shear, a and warp amplitude. For small a and small warps, 
V 2 is significantly stronger because of the resonance between disc orbits and forcing 
in a warp - the resonance is controlled by a damping, so V 2 is inversely proportional 
to a. The angular momentum transport is mainly through Reynolds stresses rather 
than viscous stresses, but this process is well described by a viscosity. 

Papaloizou & Pringle (1983) and Ogilvie (1999) assume that a acts isotropi¬ 
cally, and this is responsible for the 0:2 1/®^ relation. However, it is not obvious 

that this local isotropy holds for a viscosity driven by MHD effects, typically due 
to the magnetorotational instability (MRJ). Pringle (1992) points out that the az¬ 
imuthal shear is secular, but that the vertical shear is oscillatory. This means that 
gas parcels displaced by AR drift further and further apart, but those at the same 
radius, but displaced by Az, simply oscillate. In the simplistic MRI picture it is then 
likely that more energy is dissipated in the azimuthal viscosity, leading to a larger 
torque. However, reduced dissipation in the vertical direction may lead to an uncon¬ 
trolled resonant response and thus a larger V 2 torque. Thus simple arguments leave 
it entirely unclear how the inclusion of MHD effects should change this picture. 

If we instead assume that the MRI drives turbulence, and that the velocity field 
is uncorrelated on scales <C H (as seems probable from numerical simulations, e.g. 
Fig. 14 of Simon et al. 2012) then the action of the turbulent viscosity seems likely 


There is a long history, in both the astronomical and fluids literatures, of modelling turbulent 
fluids with effective viscosities. 
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not to care about the direction of the shear on which it acts. This supports the 
isotropic a assumption. 

The question of whether a is isotropic has been studied with numerical (Torkels- 
son et al., 2000) and analytical (Ogilvie, 2003) techniques. Both of these investiga¬ 
tions conclude that the isotropic a assumption is valid for an MRI-turbulent disc. 
There have also been attempts to interpret the dynamics of warped discs by using 
observations. These currently support the isotropic picture (King et al., 2013). 

The viscosity assumption in the diffusive and wave-like cases are effectively 
the same. As shown above, the term which describes a damping of propagating 
waves in (34) can be readily understood as exactly the V 2 torque in (27). The physics 
of both torques is the same - a damping controls the resonant response to radial 
pressure gradients which force the gas at the local orbital frequency (when the disc is 
Keplerian). When a < H/R the disc responds by propagating a wave which damps 
non-locally through a. For a > H/R this can instead be thought of as the wave 
being damped locally, and thus acting in a diffusive manner. 


3 Conclusions 

Astrophysical discs are often warped. This occurs when an initially planar disc is 
misaligned to a component of the potential (e.g. the spin of a black hole, or the 
orbit of a companion star) or when the disc becomes unstable to tilting by e.g. 
tidal (Lubow, 1992) or radiation (Pringle, 1996) effects. The main extra physics 
in a warped disc is the inclusion of a second component of the orbital shear which 
drives an oscillatory pressure gradient in the fluid. In a near-Keplerian disc this re¬ 
sults in a resonant response which communicates the tilted component of angular 
momentum. As this response is damped by an a viscosity, it depends inversely on 
that viscosity. This behaviour allows the dynamics of warped discs to split broadly 
into two regimes; one where the warp propagates diffusively and the other where it 
propagates as waves. 

The equations of motion in the diffusive case are well understood, in both the 
linear and nonlinear regimes, when the viscous stress is proportional to the local 
pressure (i.e. isotropic). The evolution equations were first derived by Papaloizou & 
Pringle (1983) and then extended to the fully nonlinear regime by Ogilvie (1999). 
This is generalised even further to include viscous dissipation and radiative transport 
by Ogilvie (2000). 

The equations of motion in the wave-like case are only really understood in the 
linear regime, again with an a damping. The nonlinear case is difficult to explore 
numerically as it requires modelling the disc with little dissipation, and thus high 
resolution. However, detailed investigation into this case is likely to be fruitful in 
revealing complex disc behaviour such as the parametric instability (Gammie et al., 
2000). This is relevant to protostellar discs where the disc turbulence may be heavily 
quenched in dead zones (Gammie, 1996). 
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There are many aspects of warped disc physics which are yet to be fully under¬ 
stood, from the inclusion of MHD, radiation warping and self-gravity to nonlinear 
hydrodynamics in strong warps. Very recently it has been shown that warped discs 
are capable of breaking into distinct planes and that this can significantly alter the 
disc’s behaviour (Nixon et al., 2012; King & Nixon, 2013). This evolution (see 
Fig. 2) has been observed in a variety of codes, with different numerical methods 
and different physics (Larwood et al., 1996; Larwood & Papaloizou, 1997; Lodato 
& Price, 2010; Fragner & Nelson, 2010; Nixon & King, 2012; Nixon et al., 2012, 
2013), so appears to be a generic feature of warped disc behaviour. It has yet to be 
shown that an MRI turbulent MHD disc is capable of breaking, but there has not yet 
been any simulation of a setup in which this would have been the likely result. 



Fig. 2 Column density projection of the tearing disc structure in a simulation following the method 
of Nixon et al. (2012) from Nixon & Salvesen (2014). 


Warped discs are relevant to a variety of astrophysical systems, including proto¬ 
stars, X-ray binaries, tidal disruption events, AGN and others. Each of these systems 
have different properties and environments which allow for a wide variety in the 
disc’s response to a warp. Warps are often induced when a disc forms misaligned to 
symmetry axes of the potential, for example the spin of a black hole or the orbit of a 
binary. However, even initially planar discs can become warped through a variety of 
effects, e.g. radiation warping, tidal torques, winds or magnetic fields. Each of these 
problems and their implications are yet to be fully understood, but with the growing 
computational power and complexity of physics employed in numerical codes it is 
likely that many questions will be answered in the next few years. 
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